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Abstract. The density of the transmission eigenvalues of Pb nano-contacts has been estimated recently 
in mechanically controllable break-junction experiments. Motivated by these experimental analyses, here 
we study the evolution of the density of the transmission eigenvalues with the disorder strength and the 
number of channels supported by the ballistic constriction of a quantum point contact in the framework 
of the Dorokhov-Mello-Pereyra-Kumar equation. We find that the transmission density evolves rapidly 
into the density in the diffusive metallic regime as the number of channels Nc of the constriction increase. 
Therefore, the transmission density distribution for a few Nc channels comes close to the known bimodal 
density distribution in the metallic limit. This is in agreement with the experimental statistical-studies in 
Pb nano-contacts. For the two analyzed cases, we show that the experimental densities are seen to be well 
described by the corresponding theoretical results. 



PACS. 5.60.Gg, 72.10.-d, 73.63.-b 

1 Introduction 

Since first transport experiments in two-dimensional elec- 
tron gases confined in semiconductor heterostructures with 
a constriction comparable to the Fermi wave length - 
quantum point contacts- where conductance quantization 
were reported [HE], much attention has been given to 
describe theoretically the quantum electron transport in 
these kind of systems. The dimensions of those structures 
are small enough that quantum interference effects can be 
observed in transport measurements, i.e., the typical size 
of the systems is smaller than the phase-coherent length. 
More recently, scanning tunneling microscope and me- 
chanically controllable break-junction (MCBJ) techniques 
have made possible the fabrication of structures at the 
atomic scale. We refer to the review article [3] for a de- 
tailed description of the progresses on those experimental 
techniques. In particular, in MCBJ experiments a notched 
wire is elongated with high mechanical precision and sta- 
bihty; thus, during this process atomic-scale contacts, or 
quantum point contacts, are formed between two wide 
electrodes. 

From the theoretical side, several aspects of the trans- 
port in quantum point contacts have been investigated 
and an extensive literature already exists. We refer the 
reader to ^ and 0] for a review of the topic. Here we 
mention some examples of the problems addressed in the 
literature where random configuration of impurities in the 
system is an important ingredient and therefore a statis- 



tical study of the transport is relevant. For instance, the 
suppression of the conductance fluctuations due to disor- 
dered reservoirs attached to a ballistic constriction was 
studied early in the nineties [51[Bj . Statistical properties of 
the conductance, resistance, and shot noise in disordered 
wires with a constriction were studied in Ref. [6] . More re- 
cently, effects of impurities around a quantum point con- 
tact on the statistics of the transmission eigenvalues have 
been investigated [7]- In single-atom contacts, it has been 
conjectured and experimentally verified that the number 
of conducting channels is determined by the valence or- 
bitals of the constriction atom 3 , 8J . 

Lately, several aspects of the electron transport in quan- 
tum point contacts have been experimentally analyzed 
using MCBJ techniques [3jl9l[T0]. For example, statistical 
properties of the transport like the density of the transmis- 
sion eigenvalues have been estimated for Pb nano-contacts 
10, llj; in these experiments, it has been pointed out that 
the transmission density for contacts with a few open 
channels is unexpectedly similar to the known bimodal 
density distribution in the diffusive metallic regime. We 
recall that in the metallic limit a large number of chan- 
nels are assumed to contribute to the transport. 

Motivated by recent break-junctions experiments where 
some statistical properties of the transmission eigenval- 
ues have been estimated for Pb nano-contacts [TUllllj . in 
this paper we calculate the density of the transmission 
eigenvalues within the framework of the Dorokhov-Mello- 
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Pereyra-Kumar (DMPK) equation [T2lfT3] . In order to ap- 
ply DMPK ideas to those break-junctions experiments, 
we assume that the constriction, or neck, formed during 
the elongation processes of the MBCJ device is smaller 
or of the order of the mean free path. This nowadays is 
achieved experimentally; in fact, the narrowest part of a 
MBCJ can be of molecular-size scale. We will consider 
that the random character of the electronic transport is 
due to the disordered wide regions at both sides of the 
contact. Finally, we suppose that the neck and the wide 
disordered region of the MBCJ device have a step geome- 
try (as in Fig. [U upper sketch). Thus, our calculations are 
based on the the solution of the DMPK equation given in 
Refs. [MlfTSlfTB] and a mapping of the transport problem 
between constricted and unconstricted wire geometries in- 
troduced in Ref. |6J . We will show that in agreement with 
the results reported in Ref. [10], our theoretical calcula- 
tions show that the transmission density can be close to 
that one in the diffusive metallic limit, even when the 
ballistic constriction of the point contact supports a few 
channels. 

In subsection 11.11 we introduce the DMPK equation 
and its solution for a single and multiple open channels. 
In the same subsection, we present the above-mentioned 
mapping of the problem of a wire with a constriction to 
that one without such constriction. In Sec. [21 we introduce 
the functions of the transmission density that we study in 
this work. Our analysis of the density start in subsection 
I2.1l with the simplest case of a point contact with a ballistic 
constriction supporting one channel. In subsection 1 2. 2 1 we 
analyze the case when several open channels are supported 
in the ballistic constriction. Also we compare some of the 
results for the transmission density reported in Ref. [TU] 
with our calculations. Finally, we give a summary of our 
study in Section [3l 



1.1 The DMPK equation 

A scaling equation for the transmission eigenvalues {ri} 
of a quasi-one-dimensional disordered system, or quan- 
tum wire, has been developed in the past [T2l[T3] . Within 
this framework, the evolution of the distribution of the 
transmission eigenvalues p{{Ti}) [or equivalently p({Xi}), 
where Xi — [1 — Ti)/Ti] with the length L of the quantum 
wire is given by the Fokker-Planck equation 
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known as DMPK equation, where N is the number of 
channels, or transverse modes, and I is the mean free path, 
while J(A) = Ylf^j \Xi — Xj\. In Eq. ([1]), we have assumed 
a system with time- reversal symmetry (/3 = 1), although 
systems without time-reversal invariance (/3 — 2) or with 
broken spin-rotation symmetry (/3 = 4) can be considered. 
Also, in the derivation of Eq. ([T]), it was assumed that the 



width of the wire is constant and smaller than its length 
L, in order to neglect diffusion in the transverse direction. 

For one channel {N — 1), the solution of Eq. ([T]) is 
given by [Hill] 
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where yo — arccosh(l + 2A) and s = L/l is the length of 
the system L in units of the mean free path /. 

The relation between the random variables A (or r) 
and the electronic transport in the wire is given through 
the Landauer formula 
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where g is the dimensionless conductance. Thus, the sta- 
tistical properties of g are governed by that ones of the 
transmission eigenvalues. 

The solution of the DMPK equation is known also for 
the general case of multiple channels and for any sym- 
metry (/3 = 1,2,4) [I4jll8j. The general solution, unfortu- 
nately, is quite complicated, but it simplifies in the metal- 
lic (1 <C s <C N) and insulating (L >> Nl ) regimes. In 
these two opposite regimes of transport, the solution of 
the DMPK equation can be written in the following form 

p(A) = lexphiJ(A)], (4) 

where Z is a normalizing factor, and H{X) is given by 
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Following Refs. [13] and [TB], we show below that we can 
use this simplified solution, Eqs. Q and (O, in our analy- 
sis. Making the change of variables A^ — sinh^ Xi, the func- 
tions u{xi,Xj) and V{xi) in the diffusive metallic regime 
are given by 
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2s * 2 
while in the insulating regime 

Uins{xi,Xj) = I sinh^ Xi - sinh^ Xj \ - In \xf ~ x'^\, 

Vins{xi) = + ^) ^2 _ iln|a;isinh2a;i| - ^InXj (7) 
Zs 2 2 

In writing Eqs. ([6]) and ([7]) we have assumed the presence 
of time-reversal invariance in the system (f3 = 1). We note 
that Eq. ([7]) differs from Eq. ^ only in the In |a;f — and 
Inxi terms. We remark that in the metallic regime a large 
number of channels, or Xi variables, contribute to the con- 
ductance, whereas in the insulating regime, the main con- 
tribution to the conductance comes from one (the smallest 
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Fig. 1. The upper sketch represents a quantum wire with 
mean free path I and a ballistic constriction in the middle. 
The wire can support A'^ channels. As explained in the text, 
the DMPK equation for such kind of wire can be mapped onto 
a DMPK equation describing a wire (as the one represented in 
the lower sketch) without the constriction, mean free path I' , 
and number of channels Nc . 



one) Xi. On the other hand, in the insulating regime, the 
average conductance is very small and Xi » 1; in this 



limit. In \xf 



and In Xi are negligible compared to the 



other terms in the insulating case Eq. ([7]). Since these 
two terms are the only ones that differ from the metallic 
solution, we thus assume that the metallic solution can 
be used as an approximated solution to DMPK equation 
valid beyond the metallic limit and for a finite number of 
channels. 

Finally, we recall that in the diffusive metallic regime, 
the density of the transmission eigenvalues defined by the 
average (^iSir — Ti)) is known to follows the bimodal 
distribution l4l[T9l: 
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with {g) = Nl/L. 



Fig. 2. (Color on-line) Density of the transmission eigenval- 
ues as defined in Eq. (|10|) for Pb nano-contacts with conduc- 
tances (in units of /h) g — 4 (solid upward triangles), 6 (solid 
squares), 10 (solid circles), and 15 (solid downward triangles). 
The solid (black) line is the bimodal distribution given by Eq. 
((S)). The inset shows the density p(r) for the same values of the 
conductances of the main frame. The data points for pnorm(T) 
and p(r) have been obtained from Refs. [10] and [TT) . 



wire geometries, in particular the solution Eq. (jUlG]), to 
study the statistical properties of transport in quantum 
point contacts; we will be especially interested in func- 
tions of the transmission density introduced in Ref. [TU] ■ 



1.2 Correspondence between constricted and 
unconstricted wire geometries 

As mentioned in the previous subsection, the DMPK equa- 
tion was derived for a wire geometry, i.e., a quasi-one- 
dimensional system of length L and width W. On the 
other hand, however, in MCBJ techniques a notched wire 
is bended to produce a small constriction between two 
wide electrodes. This constriction can be of the order of 
the mean free path of the sample, i.e., a ballistic con- 
striction. The problem of transport through a quasi-one- 
dimensional wire with a constriction can be mapped on to 
a problem of a quantum wire without the constriction, as 
it was shown by Beenakker and Melsen [B]: it turns out 
that the DMPK equation for a wire with a ballistic con- 
striction with mean free path I and N number of channels 
(Fig. [2 upper wire) is equivalent to that one for a wire ge- 
ometry (without the constriction. Fig. [Tl lower wire) with 
mean free path I' given by 

I' = 1{N+1)/{N, + 1), (9) 

where A'';; is the number of channels determined by the 
ballistic constriction in the wire. We note that a step- 
constriction geometry is assumed in the model. This map- 
ping allows us to use the known results for unconstricted 



2 Density of the transmission eigenvalues 

The density of the transmission eigenvalues p(t) were es- 
timated in Ref. [T^ by analyzing several samples of Pb 
nano-contacts. To be specific, the functions Pnoim{T) and 
q{T) defined as 

Pnorm(T) = "^^(■^)' (^'^^ 

and 

q{r) = ^ l^r'r'pir'), (11) 

were introduced by Riquclme et. al. p!Ollllj . We find con- 
venient and motivating to reproduce here some of the re- 
sults of Refs. [inj and [TT]. In figure 2, we show Pnoim(''') 
and p{t) for Pb contacts with conductances ranging from 
4 to 15 (in units of the conductance quantum 2e^/h ) [20j . 
We can observe in Fig. [5] that yOnorm(T) comes close to the 
metalHc limit result [solid hue, Eq. ([8])], even for contacts 
with small conductances. This was the unexpected result 
pointed out by the authors of Ref. [TUj . 

Coming back to the theoretical model, according to 
the mapping introduced above and assuming that Nc is 
the number of channels determined by the constriction in 
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the wire, the transmission density is given by 

p{r)^(^f;^dir-n)y (12) 

2.1 One-channel case 

We start with the simplest case of a quantum point contact 
with one open channel (iVc = 1). It is instructive to start 
with this case, although no experimental data are available 
to compare to. Using that r = 1/(1 + A) in Eq. ((2), we 
immediately obtain for the density: 
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where j/o = arccosh(2/r — 1). We also have defined the 
parameter ^ = (A^ + l)nn Eq. ((T5)) . 

Performing numerically the integral in Eq. (|13p , we ob- 
tain the density Pnorm{T) and the cumulative function g(r) 
given by Eqs. (fTO|) and (lll|) . respectively. We note that for 
the single channel case g = t. In Fig. O we show the evo- 
lution of Pnorm(''') as the value of the disorder parameter 
L/^ is changed. As expected for L > ^, small values of t, 
or conductance g, are favored, while as the ratio L/^ is de- 
creased, the density is accumulated at values of r near 1. 
As a reference, we have included the bimodal distribution 
(solid line) in Fig. [31 The function g(r) is plotted in the in- 
set of figure[3]for the same values of L/^ in the main frame. 
From Eq. ((8]), in the metallic limit: qnict{T) = 1 — \/l — r 
(solid hue in the inset of Fig. [3]). In all the cases, Pnorm(T) 
and g(r) clearly exhibit discrepancies with respect to the 
diffusive metallic limit, in contrast to the results for the 
multichannel case, as we will see in the next subsection. 



2.2 Several-channel case 

Let us now assume that the ballistic constriction can sup- 
port several channels Nc- Therefore from the expression 
for the density [Eq. (|12p ] and using the join probability 
p{{xi}) given by Eqs. ([4]- [6]), we can write the density of 
the transmission eigenvalues as 



exp 



1<] l<i<j Kj 



(14) 



where x = arcsech-y/r, u(xi,Xj) is given in Eq. and 
{N+1) I , 1 



ln(xi sinh 2xi). 



(15) 



We will see that the density of the eigenvalues r„ given 
by Eq. approximates to the diffusive metallic regime 
[Eq. ([8])] even for a relative small number of channels Nc. 




Fig. 3. Main frame: Density pnorm('r) for A'^c = 1 at dis- 
order parameter 2L/^—2 (dashed line), 2/5 (dotted line), 
1/5 (dashed-dotted line) with average conductance (g) — 
0.26,0.71,0.83, respectively. The solid line is computed from 
Eq. (|8]) and corresponds to the diffusive metallic regime. Inset: 
The cumulative function ^(t) given by Eq. (|lip for the same 
values of the parameters 2L/^ as in the main frame. 
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Fig. 4. (Color on-line) Density pnorm(r) for Nc = 3,4, and 
6 (dotted red, dashed green, dashed-dotted blue lines, respec- 
tively). In the inset g(r) is shown for the same values of N^ 
in the main frame. For all cases, the average transmission (r) 
is approximately ((r) = 0.49), while the average conductance 
(3) = 1.5,2,2.8 for A^^ = 3,4, and 6, respectively. The solid 
(black) line is the bimodal distribution in the metallic limit, 
Eq. dH). 



We next present the results for Pnorm(''') and g(t) [Eqs. 
(fTO)) and pT|) ] obtained from the above expression for p(t) 
[Eq. (jl4p ]. where the integrals were performed numerically. 

Let us start showing the evolution of the density of 
eigenvalues with the number of channels A^c • In figure [H 
Pnorm(''') and q(i) are plotted for N^ = 3,4, and 6; these 
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Fig. 5. (Color-online) Evolution of the density pnorm(T) with 
the disorder parameter I//^. In the inset, g(r) is plotted for 
the same values of L/(_ in the main frame. In all the cases the 
number of channels is fixed to Nc = 5. The solid (black) line 
is the bimodal density in the diffusive metallic regime. 



cases are chosen to have the same average value (r). We 
can observe that as the number of channels is increased, 
Pnorm(''') concentratcs at small values of r and near 1, 
i.e., the density evolves to the bimodal distribution plot- 
ted in solid line in the same Fig. 21 In other words, as 
the conductance of the ballistic constriction increases, the 
distribution of the transmission eigenvalues approximates 
quickly the bimodal distribution. On the other hand, the 
cumulative function ^(t) (inset in Fig. 3]) also comes close 
to the result in the metallic limit gmet(T) (solid line). In 
fact, we note that the different curves for q{T) might look 
closer to the metallic case (jmct (t) than the corresponding 
cases for /9norm(''') to Pmct(T) (the bimodal distribution). 
This is due to loss of details of the density Pnorm {t) when 
this function is integrated [Eq. pTjl ] to obtain ^(t). 

Now we show the evolution of the density with the 
ratio L/£^. In figure [5l we plot the results for Pnorm{T) 
and ^(t) (inset) at three different values of the disorder 
parameter L/^ for a fixed value of the number of chan- 
nels Nc = 5. For the smallest value of L/^ (dashed line), 
PnormlT) shows a peak at r's near 1, like the bimodal dis- 
tribution (solid line); in fact, the density PnormiT) looks 
similar to this metallic limit case, except at small values 
of T. As the ratio L/£^ is increased, i.e., going to the insu- 
lating regime, the peak at small values of t grows rapidly 
and /CnormCT") comcs close to the bimodal distribution at 
small values of t. An increase in peak height at small r's, 
however, causes a detriment to the peak at r near 1, as 
one might expect: roughly speaking, we can say that the 
evolution of yOnormlT) with the ratio L/^ described above is 
consequence of the normalization of the density function. 
Regardless of this behavior, for relatively small values of 
the ratio L/£^, the main differences of Pnorm(''') with re- 
spect to the diffusive metallic limit are seen only at small 
values of r, despite the small number of channels of the 
constriction. 
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Fig. 6. (Color on-line) Comparison between the experimen- 
tal results (solid symbols) for the density p(r) for Pb nano- 
contacts and the theoretical calculations (dashed and dashed- 
dotted lines). The conductances (in units of e^/h) of the ex- 
perimental data points are g = 6 (solid squares) and 10 (solid 
circles), also shown in the inset Fig. [2] The theoretical results 
are chosen to have average values (r) similar to the experimen- 
tal ones, which we have estimated from the experimental data 
in the region (0.14 < r < 1), Fig. [21 (r) r; 0.4 (solid squares), 
0.5 (solid circles). The theoretical curves are normalized to the 
area of the experimental density in the region (0.14 < r < 1). 
A good agreement is seen between the theoretical densities and 
experimental data points in the region where experimental re- 
sults are reported. Inset: density pnorm(''") for the same cases in 
the main frame. The solid (black) line is the bimodal density 
in the diffusive metallic limit. 



Finally, we compare some of the experimental results of 
Refs. [lOj and [TT] (Fig.jJ) with theoretical predictions. We 
point out that the theoretical model depends on the num- 
ber of channels Nc and the disorder parameter L which 
we do not know from the experiments of Ref. [lOj. Thus, in 
order to make the comparison, we have extracted the av- 
erage transmission eigenvalue (r) from p(r) plotted in the 
inset of Fig. |2l For the couple of chosen cases, solid squares 
and solid circles in Fig.j^l we have estimated: (r) « 0.4 and 
0.5, respectively. We also note that these averages values of 
T are calculated in the interval of r (0.14 < t < 1 ), region 
where experimental data points have been reported. From 
our calculations, the cases Nc = 6 and 7 with L/^ = 2/7 
and 1/8, respectively, give mean values (r) similar to the 
experimental ones. In Fig. |S1 (main frame) we compare the 
results for p{t): the dashed line corresponds to A^c = 6, 
while the dashed-dotted line is the result for Nc — 7. In 
order to compare with the experimental results, we also 
have normalized the theoretical density to the area of the 
experimental density in the interval (0.14 < r < 1). We 
can see that the experimental density is well described by 
theory. In the inset of figure |6l we also compare the nor- 
malized density Pnorm(''") for the same cases in the main 
frame. The solid line in the inset corresponds to the bi- 
modal distribution. 
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3 Summary 

We have studied the evolution of the transmission eigen- 
value density of quantum point contacts with the disor- 
der strength and the number of channels of the ballistic 
constriction. Our analysis of the transmission density is 
based on an approximation P^SlfTB] to the general solu- 
tion of the DMPK equation and a mapping of the prob- 
lem of transport through a quantum wire with a ballistic 
constriction to the problem of transport through a wire 
without the constriction [6] i.e., the latter problem is de- 
scribed by the known solution of the DMPK equation. We 
have shown that the density of the transmission eigenval- 
ues comes close to the known bimodal distribution in the 
diffusive metallic regime even for a small number of chan- 
nels supported by the constriction of the quantum point 
contact. This is in agreement with the unexpected result 
pointed out in Ref. '10], where some statistical properties 
of the transmission eigenvalues for Pb nano-contacts were 
analyzed by using MCBJ techniques. In order to apply 
the DMPK approach to those experiments, we consider 
that the neck formed during the elongation process in a 
MCBJ device is a ballistic constriction which is in between 
two wide disordered regions. We do not have experimental 
information of each parameter that enters in the theoret- 
ical model, but by comparing theoretical and experimen- 
tal results with similar average transmission values (r), 
we have seen that the analyzed experimental densities are 
well described by the corresponding theoretical densities. 
The trend of the transmission density seen in Pb contacts 
experiments as the number of channels is changed is also 
in agreement with the theoretical results presented here 
when the number of channels of the ballistic constriction 
is varied. 
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